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New features of modulational instability of partially coherent light; 
importance of the incoherence spectrum 
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It is shown that the properties of the modulational instability of partially coherent waves propagating in a non- 
linear Kerr medium depend crucially on the profile of the incoherent field spectrum. Under certain conditions, 
the incoherence may even enhance, rather than suppress, the instability. In particular, it is found that the range 
of modulationally unstable wave numbers does not necessarily decrease monotonously with increasing degree 
of incoherence and that the modulational instability may still exist even when long wavelength perturbations are 
stable. 
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The modulational instability (MI) of coherent constant am- 
plitude waves in nonlinear Kerr media is one of the most 
fundamental phenomena resulting from the interplay between 
nonlinear phase modulation and linear dispersion/diffraction 
and has attracted much interest over many years Q,|3>|3|. Re- 
cent advances in the area of nonlinear optics, in particular new 
results concerning the nonlinear propagation of partially inco- 
herent light and the advent of incoherent soli tons 0], have 
prompted a revisal of this issue during the past decade. The 
general understanding that emerged from these studies is that 
the wave intensity threshold for the onset of the MI is in- 
creased by the wave incoherence. With this picture in mind, 
it is remarkable that in a recent investigation of the transverse 
instability (TI) of solitons, (sll, it was found that the range of 
modulationally unstable wave numbers did not monotonously 
decrease with increasing degree of incoherence. In fact, it first 
increased until eventually it started to decrease. 
Inspired by this result, we consider, in the present work, the 
problem of the modulational instability of partially coherent 
waves in more detail and show that the picture is more com- 
plicated than previously thought. In order to simplify the 
analysis and to bring out clearly the new features, the anal- 
ysis is carried out for the longitudinal modulational instabil- 
ity. We find that the effect of the incoherence on the MI is 
sensitive to the profile of the incoherent power spectrum. For 
the often used assumption of a Lorentzian incoherence spec- 
trum, the range of unstable wave numbers does indeed de- 
crease monotonously with increasing degree of incoherence, 
whereas e.g. for a Gaussian spectrum, the range first increases 
and then starts to decrease monotonously. This result agrees 
well with the unexpected feature observed in Q] . Also, several 
other subtle effects are shown to be possible. In particular, it is 
found that the threshold for the MI to be completely quenched 
is not necessarily associated with the long wavelength limit. 
Modulations may be stable in this limit, but still be unstable 
for finite wave numbers. This implies that the threshold for 
total quench can not, in a general case, be determined by sim- 
plifying the analysis to considering only the long wavelength 
limit as is done in e. g. 0,0] ■ 

The starting point of our analysis is the normalized nonlinear 
Schrodinger equation describing the one dimensional prop- 



agation of a partially incoherent wave in a dispersive (or 
diffractive) nonlinear medium, viz 



.dtp Id'^i' 



(1) 



where the bracket, (. . .), denotes statistical average H. This 
equation is valid under the assumption that the medium re- 
sponse time is much larger than the characteristic time of the 
stochastic intensity fluctuations. 

The modulational instability of small perturbations of the 
corresponding steady state solution has been analyzed using 
different, but equivalent mathematical formalisms. An 
analysis based on the formalism of the correlation function, 
iIsIFtI IicIi or on the Wigner approach 01, results in the disper- 
sion relation 



po{p-k/2)-po{p + k/2) 
kp + Z7 



dp = 1, 



(2) 



with pq being the Wigner distribution function of the un- 
perturbed cw wave. However, using the transformations: 
p + k/2 ^ 9,p — k/2 = 6', Eq.(|2|i can be expressed as 



+ 00 



Po{o')de' 



, k{e' + k/2 + ij/k) 

This is then easily rewritten as 



+ 00 



Po{9) de 



h{9 - k/2 + ij/k) 



Po{o) de 



fc4/4+ (a-6i + 7)2 



= 1 
(3) 



(4) 



which is the same expression as the Coherent Density ap- 
proach Ol, provided we identify ^0(6*) = G{6) with 
being the averaged normalized field intensity of the station- 
ary state and G{9) being its normalized angular spectrum. 
Throughout this paper we will use both expressions for the 
dispersion relation interchangeably, since some parts of the 
analysis are most conveniently handled by one approach, and 
some by the other 

An explicit analytical solution of the dispersion relation, 
Eq.Q or Eq.0, is possible only for some particular inco- 
herence spectra. Specifically, in the case of the Lorentzian 
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spectrum G{9) = 9q/[k{0'^ + ^olL where 6q is the width of 
the spectrum, it has been shown jsl], that the solution can be 
expressed as 



7(fc,^o)=7o(fc)-fceo 



(5) 



where 70 (fc) is the growth rate in the coherent case {G{9) = 
(5(6')) i. e. 



7o(fc) = k^A^ - fc74, 



(6) 



with > This analytical result shows explicitly that 

the effect of the incoherence, provided it is large enough, is to 
suppress the MI for any value of the perturbation wave num- 
ber However, as will be demonstrated in this letter, this result 
depends crucially on the form of the incoherence spectrum 
and is not a general feature of the MI of partially incoherent 
light. 

The restricted generality of the result expressed by Eq.ljSj can 
be directly inferred by studying in more detail the proper- 
ties of the cut-off wave number, kc, i. e. the value of k at 
which the growth rate vanishes. According to Eq.(|5jl, valid 
for the Lorentzian spectrum, kc is shifted monotonously to 
the left (decreased) with increasing degree of incoherence, Oq: 
k^ = 4(A^ — 9q). For the case of a general coherence spec- 
trum, kc is determined by the following equation: 



A' 



+ 00 



G{9)d9 



1 



(7) 



where p.v. denotes the principal value. When the power spec- 
trum is well localized, in the sense that its rms-width is much 
smaller than the cut off wavelength, the contributions from the 
zeros of the denominator are negligible and the denominator 
can be expanded to yield 



1 



kl/A 



1 



kl/A 



d9. 



(8) 



Keeping only the first two terms of the expansion one obtains 
an approximate solution for the effect of partial coherence on 
the cut-off wave number as follows: 



(9) 



where {9'^) = j 9^G{9)d9. Thus kc is found to increase for 
increasing degree of incoherence, which at first sight seems 
to be in contradiction to the behavior of the previously found 
exact solution for the Lorentzian spectrum, cf Eq.(|5|l. How- 
ever, the Lorentz spectrum is not well localized in the sense 
defined above, since the value {9'^) does not exist. A first in- 
dication of such an incoherence-induced increase of the range 
of modulationally unstable wave numbers was observed by 
Torres et al. jsl, in a numerical study of the transverse in- 
stability (TI) of soliton structures using a Gaussian spectral 
distribution. It should be emphasized though, that the effect 
will occur for transverse as well as for longitudinal modula- 
tional instabilities. Although the Lorentz spectrum proper is 



not well localized, it can easily be made so by considering the 
bounded Lorentz spectrum i. e. 



G{9, 9o, 9m) 



1 



cw{e^-\9\) (10) 



where W{x) = if a; < 0, W{x) = 1 if x > 0, 6*™ is the 
boundary of the spectrum and C is a normalization constant. 
For this spectrum one can explicitly show that, depending on 
the value of Om, the cut-off shift may either increase or de- 
crease with increasing 9q. 

Another perturbative solution of the dispersion equation for 
the MI is possible in the long wave limit when k <C ^^o. where 
6'o is the characteristic width of the spectrum 0(6*). In this case 
one can use a Taylor expansion around p of the numerator in 
Eq.Q. Introducing F = 7/fc we obtain 
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-A' 
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G'{P) 
P + iT 



dp- 
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24' 



-A^k" 



+ 00 



G"'{p) 
ID + iF 



dp. (11) 



This simplified approximation (however without the last term 
on the right hand side of Eq. Jl It ) was used previously |7| to 
analyze the threshold condition for the suppression of the MI. 
Actually, assuming G{p) = G(-p) and F = F* < 6*0, the 
first integral can be approximated as 



G'{p) 
P + iT 



dp = — Ji + ttFDi, 



where 



Ji 



+ 00 



G' 



dp, 



Di 



P 



(12) 



(13) 



p=0 



Within this approximation the dispersion relation becomes 



ttDiF = Ji- a- 



(14) 



The equality Ji ~ A^^ can thus be taken as determining the 
threshold for MI development. While this is true in the limit of 
vanishing k, it is clear that by taking into account also the next 
term in the expansion Eg. dl II . instead of Ea.( ll4t one obtains 
the solution: 



ttDiT = ,h- A-^ - ^k^J3, 



where 



+00 (~im 



V 



■dp. 



(15) 



(16) 



This analysis is valid provided D\ 7^ 0, in other case a full so- 
lution of Eg. dl \\ is needed. Performing the calculations one 
finds that for a Gaussian or a Lorentzian incoherence spec- 
trum, the values of D\, J\, and J3 are all positive quantities. 
Thus in such cases, Eq. (ll4t provides a sufficient condition for 
the suppression of the instability, as found in the works of 
Anastassiou et al. 0] and Soljacic et al. 0]. Nevertheless, 
there may exist spectra, for which the factor J3 is negative. 



3 



implying there is positive growth of the perturbation for finite 
k, despite the fact that solution of the EQ.( ll4t (the long wave 
limit) gives < 0. Consequently, the definition for the sup- 
pression of the modulational instability as the threshold value 
given by the long wavelength limit approximation is not ap- 
propriate. 

The simplest illustration of the ambiguity of the threshold con- 
dition based on the long wavelength limit can be given by an- 
alyzing a rectangular spectrum profile with G{p) = l/26',„ 
for p in the interval —9,n<P< djn, whereas G{p) = 
outside this interval. For this simple spectrum, the integral in 
Eq.@ may be evaluated exactly and the following dispersion 
relation is obtained 



ta.nh{ke„JA^) 



(17) 



In the limit k6„ 
imated to read 



<C A^, the dispersion relation can be approx- 



3^2 



(18) 



The shift of the cut-off wave number kc 
mately by 

kl^A{A'+0l/3), 



IS given approxi- 



(19) 



provided dm ^ A, which is in full agreement with Eq.Q. 
However, the rectangular spectrum exhibits one more unex- 
pected and very important feature: It is evident that even if 
7^ < in the limit as k 0, 7^ may become positive for 
finite k. In particular, when k ~ 2dm, the growth rate of the 
perturbation is positive, 7^ > 0, independently of the spec- 
trum width 9m, in fact 



7 



16 9t 



exp(4 02^M2)_i 



> 0. 



(20) 



Therefore, even a very high degree of incoherence does not 
completely suppress the modulational instability when the 
spectrum is rectangular, cf. FigQ. Even when 9m ^ A 
and the MI is strongly suppressed within the long wave range, 
there still exists a "resonance" region of instability for wave 
numbers k around 29m given by : 



\k-29m\ <46',„exp 



A^ 



(21) 



It is interesting to note that a similar phenomenon was found 
in fT3l . in an investigation of the effect of a nonlocal nonlin- 
earity on the MI in a focusing Kerr medium. In particular, it 
was shown for a rectangular response function, that the insta- 
bility growth rate first decreased for increasing width of the 
response function, but ultimately for large widths, new modu- 
lationally unstable bands appeared at finite wave numbers. 
Since a rectangular profile is a very ideal and a bit artifi- 
cial form of the spectrum, we consider another example of a 
well localized spectrum in the form of a modified Lorentzian, 




FIG. 1 : The effect of increasing rectangular spectrum width 9m on 
the MI. The parameter 9m runs from 6m ~ (the top most curve) in 
increments of 0.25 to 9m ~ 3 (the bottom most curve). As can be 
seen, the instability is never completely suppressed. 



which exhibits similar properties to those of a rectangular pro- 
file. 



G{9) 



V2 91 



(22) 



Even in this case the integral of Eq.@ can be calculated in 
explicit form to yield 



+00 



G{9)d9 



(7 + iA:0)2 + (fc2/2)2 
W4-^ + (7+^) +2(7+^)^ 



(23) 



^2 / \/2 



Using this expression in the cut-off condition, 7^ = 0, we 
obtain the following result for the cut-off wave number: 



kJ.^2A^<l± 



\ 



(24) 



When 1 < (6*0/^)2 < (1 + V2)/2, the equation has two pos- 
itive roots implying that the MI is not completely suppressed 
provided this condition is fulfilled. On the other hand, 72 > 
only within a limited range of wave numbers, viz 



fc2 
2^2 



< 



\ 



A 



(25) 



In order to further illustrate the subtlety of the interplay be- 
tween the partial incoherence and the instability drive, we 
consider a multi-carrier case where the field consists of many 
mutually incoherent, but individually partially coherent waves 
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with a spectrum given by G{9) = 
persion relation then becomes 



Gnid-en)de 



Or,)- The dis- 



(26) 



For simpHcity we consider the particular case of equally sep- 
arated carriers, i. e. 9n = an where n is an integer, with each 
carrier having a Lorentzian phase spectrum of the same width 



din) 00 



" Ti\{B - any + el] ■ 

The dispersion relation, Eq.(l26>. can then be written as 



(r + ifcan)2 + fc4/4 



1, 



(27) 



(28) 



where F = 7 + kO^- For a symmetric spectrum, i.e. when 
g{k) ~ g{—k), this dispersion relation can be rewritten in 
terms of real functions: 



X + {a + na){a — na) 
+ (a + na)2] [X + {a- naf 



1 (29) 



where X = V'^/k'^ 



k/2 and the coefficients g{n) are 



normalized to unity, g{n) = 1. When X has a small pos- 
itive value, the sum in Eq.(l29> has multiple resonant values at 
a,i = an, n ^ 0. Consequently, one should expect the ex- 
istence of small positive roots X of Eq.(l29> in the vicinity of 
those resonances, independently of the structure of the spec- 
trum envelope g{n). That this indeed is the case is illustrated 
in Fig. (13, which shows the result of a numerical evaluation 
of the sum for particular parameter values. On the other hand. 




0.1 r 




FIG. 2: Numerical solution of EQ. <29t for multi carrier operation. 
The number of carriers is 10, a — 0.2, and no = 0.8. Top: Main 
growth curve together with the small instability islands stemming 
from the separate carriers. Bottom: A blow up of the instability 
curve. Notice the resonant peaks occurring beyond the region of the 
main curve. 



outside of these resonant bands, the summation can be trans- 
formed into an integration over n provided the spectrum g{n) 
is dense enough i. e. a <C A \g{n + 1) — g{n) \ <C g{n); 



+00 



y_oo ■'^^"^ [X + {a + nay][X + {a - na)^] 

(30) 

This equation coincides with Eq.@ if 17(71) dn is replaced by 
G{9) dO and na is changed to 9. It corresponds to the MI of 
a continuous spectrum which coincides with the envelope of 
the actual spectrum. For example, when 



9{n) = - 



no 



1 



Trn^ -I- riQ ' 



(31) 



where no ^ 1, Eq.dSOt is reduced to the following well 
known expression: 



(r + kanof = k^{A^ - k^/A). 



(32) 



Thus, within this approximation, the MI is suppressed for all 
wave numbers provided ano > A. However, as shown above, 
this is correct only outside the resonance bands, i.e. the lo- 
calized regions around each 2an (n ^ 0), where F^ > 0, 
independently of the width of the envelope (no). 
To conclude, we have investigated the role of the incoherent 
spectrum profile on the properties of the modulational insta- 
bility. The gain curve of the instability may smoothly shrink 
in amplitude and cut-off wave number with increasing degree 
of incoherence, as is the case for a Lorentzian profile. How- 
ever, it may also initially expand into the wave numbers which 
are stable in the coherent regime as is the case for a Gaus- 
sian profile of the spectrum. When the spectrum is rectangu- 
lar, we have shown that under the special resonance condition 
k ~ 29m, the MI cannot be suppressed completely no mat- 
ter how strong the degree of incoherence. Using the modified 
Lorentzian profile we demonstrate that the long wavelength 
threshold definition at which the MI is suppressed due to par- 
tial incoherence is in fact profile dependent. For this case, is- 
lands of positive growth rate emerge at higher wave numbers. 
Lastly, our analysis of the modulational instability in the case 
of multi-carrier operation demonstrates new additional struc- 
ture in the gain curve, where besides the main lobe, there also 
exist smaller peaks surrounding the discrete carrier phases. 
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